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Abstract 

A crossover from d to d—1, and then back to (/-dimensional critical behavior 

is argued to be a generic feature characterizing ordering in a d-dimensional 

superlattice composed of atomically thick films of two ferromagnets. The 

crossover leads to anomalous changes in the amplitudes of critical singularities. 

In d = 3 Heisenberg and XY superlattices large scale critical fluctuations 

persist over a wide temperature range. 
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I. INTRODUCTION AND THE STATEMENT OF THE PROBLEM 



Recent progress in creating artificial magnetic heterostructures ]IJ gives an impetus to 
the theory of phase transitions in multilayers. While the related experimental and theoretical 
efforts have been so far mostly concentrated at microscopic effects, this paper addresses the 
problem of universal features associated with large scale ordering in superlattices. Suppose 
the superlattice is built of atomically thick layers of two magnets, which taken separately 
both order through second order phase transitions at two respective bulk critical points 
T cl > T c2 . Suppose also for simplicity that both components are ferro- magnets and that 
the critical behavior at their bulk Curie points is similar, i. e. both components belong to 
the same universality class. The question then is: at what temperature T c will long range 
ferromagnetic order appear in the superlattice and what critical behavior will be observed at 
that temperature? One can reasonably expect that in the limit of thick layers the answers 
to these questions will be largely universal, depending only on the bulk critical properties 
of the two components and the geometry of the superlattice; these universal aspects of the 
problem represent the subject of this paper. 

We will consider the simplest superlattice geometry (Fig. 1) constructed of two elemen- 
tary building blocks: slabs, or layers, of two ferromagnets, 1 and 2, of finite thickness Li and 
L 2 , respectively. The slabs are stacked periodically in the z-direction, so that L = L\ + L 2 is 
the period. The system is homogeneous in the remaining d! = d — 1 dimensions. While the 
dimensionality d = 3 (Fig. la) is naturally the most interesting one in view of experimental 
applications, the planar, d = 2, superlattice geometry can be quite conceivably realized by 
cutting a thin film out of a three-dimensional multilayered sample (Fig. lb) or by deposition 
of a magnetic film on a substrate, which is itself cut out of a superlattice. Extension of our 
consideration to d = 2 is important because of analytical tractability of two-dimensional 
models. Being interested in the long wave length aspects of the problem, we will consider 
the limit of atomically thick slabs: L\ ^> a\, L2 > a 2 , where a\^ 2 are the thicknesses of 
elementary, molecular layers of the two components, respectively; correspondingly, the in- 
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terfaces separating two subsequents slabs do not have to be atomically smooth. Naturally, 
universality extends applicability of the results obtained below beyond this particular mag- 
netic model; among other interesting applications are a superfluid in a periodic matrix and 
a superfluid film deposited on a periodic substrate. 

Universal critical phenomena related to three different geometrical elements present in 
a superlattice have been extensively studied theoretically and, to a variable extent, ex- 
perimentally. First, at a single interface separating half-infinite near-critical, say \t\\ = 
\T — T c i\/T c \ <C 1 and disordered, \t 2 \ = \T — T c2 |/T c2 = 0(1), magnets, one expects ordi- 
nary surface critical phenomena to take place at the near-critical side of the interface. It has 
been established theoretically 0,^] that the critical exponents and amplitudes characteriz- 
ing surface values of various observables are quite different from those in the bulk, as well 
as from the a 1 '- dimensional ones. However, experimental confirmation of those theoretical 
findings has turned out to be difficult. Note that near T c2 , where the other half of the single 
interface model goes through the criticality, the first one is already ordered and imposes a 
magnetic field acting at the surface of the second component. Therefore, one expects normal 
surface critical behavior ||,|5| to take place at the interface near T c2 . 

Second, near the bulk critical temperature T CJ - of the j-th, j = 1,2, component a thick, 
Lj ^> Oj, single slab of that component exhibits a dimensional crossover from the critical 
behavior characteristic of bulk, d-dimensional, samples to the d' = d — 1-dimensional critical 
behavior, as observed in thin films of the same material 0. We will refer to this crossover 
as bulk-to-film below; it has been observed experimentally |7j. Finally, in a system of thin 
layers connected to each other by very weak interlayer bonds, near the critical temperature 
T' c of a single layer one expects a crossover from d'- dimensional to the bulk, <i-dimensional 
behavior, which we will call film-to-bulk crossover below. This type of dimensional crossover 
has been well experimentally documented in a variety of magnetic H and non-magnetic 
layered systems |5). 

One source of motivation for studying critical phenomena in multilayers is that the 
significantly increased surface-to-volume ratio may improve chances of observing features 
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characteristic of the d = 2-dimensional and, especially, the elusive d = 3-surface critical 
behavior with respect to experiments involving single thin film or semi-bulk samples. How- 
ever, finite separation between interfaces in a superlattice leads to correlations between 
fluctuations in different layers and at different interfaces. One therefore expects surface 
and two-dimensional scaling to be limited to a certain range of length scales, while even- 
tually cooperative phenomena involving fluctuations on scales larger than the period of the 
superlattice will dominate the criticality. A crucial insight comes from the exact solution 
available for a ferromagnetic d = 2-Ising superlattice, as modelled by a planar Ising model 
composed of alternating strips of two components flQfl . Recent analysis |nj demonstrated 
that the anomalous changes in the critical amplitudes and other interesting features of that 
model can be simply explained by reentrant dimensional crossover, as the temperature ap- 
proaches the critical temperature of the superlattice the long wave length properties of the 
system are successively dominated by fluctuations characteristic of the the d = 2-Ising, then 
d = 1-Ising, and then again d = 2-Ising critical behavior. 

Here we generalize and explore this picture by constructing and analyzing a renormal- 
ization group (RG) flow which provides a scaling description of criticality in experimen- 
tally interesting three dimensional superlattices as well as in the Heizenberg (isotropic) and 
XY (easy plane) universality classes. In the proposed scenario as the scaled temperature 
t — (T — T c )/T c , or ti — (T — T cl )/T cl , decreases, the critical behavior of the superlattice 
displays the features of the three simpler systems listed above (Fig. 2): first, as T cl is 
approached, the original superlattice ((i) in Fig. 2) can be approximated by sequence of 
thick layers of the first component weakly coupled through the 2-layers ((ii) in Fig. 2). At 
larger values of t this coupling can be neglected and the layers display regular bulk critical 
behavior with the ordinary surface critical behavior at the interfaces. Then, as the cor- 
relation length in the first component becomes comparable to the thickness L 1; each slab 
of the first component exhibits direct bulk-to-film, d — > cf-dimensional crossover. At still 
smaller values of t, the thickness of these slabs is irrelevant and the system becomes equiva- 
lent to the third paradigm mentioned above: a weakly coupled layered system ((iii) in Fig. 
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3). The reverse, film-to-bulk crossover finally takes the system back to the uniform, bulk 
d-dimensional behavior, characterized however by high degree of anisotropy ((iv) in Fig. 2). 
This reentrant crossover behavior leads to dramatic changes in the critical amplitudes, as 
shown in Fig. 3 for the d = 3-Ising class: the critical exponent characterizing divergence of 
a physical quantity, for which the specific heat per unit volume C has been chosen in Fig. 3, 
changes two times, as t goes to zero, between the values a and a' characteristic of the d = 3 
and d' = 2-dimensional critical behavior, respectively. The relative widths td (subscript 
D denotes direct, bulk-to-film crossover here) and tr (subscript R stands for the reverse, 
film-to-bulk crossover) of the temperature domains where the planar Ising and the reentrant 
bulk critical behaviors are observed will be calculated below; they provide convenient scaling 
combinations in terms of which a proper description of criticality in superlattices is achieved. 

Two further points have to be added to this preview of the paper. First, a crucial element 
of our qualitative picture is that different layers of the first component remain effectively 
decoupled while the bulk-to-film crossover (from (ii) to (iii) in Fig. 2) happens within them. 
This turns out to be guarantied by the fact that the subsequent 1-slabs are coupled through 
the surface spins whose correlations are much weaker than of those in the bulk. Thus the 
surface critical behavior, which is not directly seen in temperature dependencies of bulk 
quantities, like that shown in Fig. 3, appears crucial to the nature of criticality in the 
superlattice; we will show that the values of surface critical exponents can thus be extracted 
from analysis of the dependence of the reentrant width tr on the thickness of the 2-slabs, 

Second, in the case of the d = 3 Ising superlattice the d! = d — 1-dimensional model, 
describing a single film of the first component, has a regular finite temperature phase transi- 
tion. Consequently, in the thick layer limit the critical temperature T c of the superlattice is 
only slightly depressed from T cl and the critical fluctuations are essentially localized within 
the 1-layers. However, in all other cases, i.e. d = 3 Heizenberg and XF-models and all mod- 
els in d = 2, the corresponding cf-dimensional universality classes lack true ferromagnetic 
long range order at finite temperatures. The interlayer coupling through the layers of the 
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weaker, second component thus are crucial to existence of long range order. As a result, the 
critical temperature in these models is strongly shifted towards T c2 and the critical domain, 
in which large scale fluctuations are active, is extended to the temperature range of the 
order of T cl — T c2 . 

Finally, a few remarks about the method of this paper. We focus here on the n-component 
vector spin systems with short range interactions. As usual, Ising, XY, and Heizenberg 
universality classes refer to the n = 1, or easy axis, n = 2, or easy plane, and n = 3, or 
isotropic, magnets respectively. The most important omission of this model, apart from 
restriction to the simplest type of the order parameter, is the absence of long-range dipolar 
forces; it seems however natural to make the first step within the simpler realm of local 
models. The actual method employed here is to combine the well known linear RG flows at 
the fixed points representing the three paradigms discussed above (Fig. 2), into the simplest 
global RG flow consistent with this linear behavior (Fig. 4). The resulting scaling theory 
then passes the test of comparison to the exact results available for the planar Ising model. 

The outline of the paper is the following: in the next section we start by taking the 
conceptually simplest case of the three-dimensional Ising superlattice. The renormalization 
group flow will be proposed and scaling forms derived. In section III we consider the planar 
Ising case which requires extension of the previous consideration to the case of the zero- 
temperature critical point describing the d! — 1 dimensional Ising model. The results will 



then be shown to agree with the exact forms available for this case [|TTJ . The extension of the 
method developed in section III will be used in the following section, where the Heizenberg 
and the XY universality classes are considered. The results are summarized in section V. 

II. THREE-DIMENSIONAL ISING SUPERLATTICE 

We start now by focusing on the d = 3-Ising universality class. Its most important 
feature is that a single thick slab of the first component is fully capable of ordering at a 
finite temperature T' cl , which continuously approaches the bulk transition temperature T cl , 
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as the thickness L\ — > oo. We thus expect that the criticality in a three-dimensional Ising 
superlattice happens close to T cl and that the relevant fluctuations are localized within the 
1-slabs. Much of what will be said below will apply however to other systems, so we will 
keep using general d and d! — d — 1, as well as the general standard notation for the critical 
exponents in this section. 

The renormalization group procedure which we use here has been extensively discussed in 
|| for the half-infinite models bounded by one free surface, and applied to the superlattice 
geometry in [jlXfl . It is based on subsequent reduction, A' — > A'e~ l , of the upper cut-off 



A' imposed on the cf-dimensional wave numbers k' characterizing the spatial variation of 
fluctuations in the directions parallel to the layers. Since the superlattice geometry is uniform 
in this plane, expansion in plane waves exp(ik'x') can be used as a basis for a perturbative 
RG of the type discussed in P,|12|. The details of the procedure will not be important to 



us here: all we need is the existence of an exact RG transformation; then according to the 
general principles [|T2| , |13| the linearized RG flow in the vicinity of the fixed points should be 
independent of the specific implementation. Note that there is no need to do anything about 
the cut-off in the ^-direction (in fact one does not need such cut-off at all): elimination of 
fast modes in the directions parallel to the layers automatically induces coarse-graining in 
the z-direction H3|, [TT1 , [T4 . 



We will now follow the RG flow probing fluctuations on larger and larger length scales 
parallel to the layers. The two crucial length scales, to which the current (running) length 
scale has to be compared, are the bulk correlation lengths £i(T) and ^{T), which can be 
(at least in principle) measured at the given temperature T in independent bulk samples 
of the two components. Near T c \ fluctuations in the second component are confined to 
scales smaller than the bulk correlation length £%(T) ~ £,2{T c i), which remains finite while £1 
diverges. Decreasing the upper momentum cutoff to £j essentially eliminates fluctuations 
in the 2-layers. On larger scales the 2-layers are adequately described by the Ornstein- 
Zernike spin density functional corresponding to the Gaussian RG fixed point describing the 
low-temperature phase of the second component. As usual, Gaussian degrees of freedom can 



be integrated out, the integration being equivalent to minimization of the Ornstein-Zernike 
functional. Such minimization is performed in the Appendix. It results in an effective 



interaction [13 1 



■Joe- L ^ j ' s T n {x)s* +1 {x)d d - l x> (1) 

between the top surface spins s^(x) of the n-th layer and the bottom surface spins s^ +1 (x) 
of the n + 1-th layer of the first component. At this point the thickness L 2 gets absorbed 
into the bare coupling constant Ji = Jo exp(— Z^/^)- A major simplification has occurred 
after this initial crossover: the RG flow has taken the original superlattice system to the 
first fixed point, FP Q , equivalent to a sample of the first component containing a periodic 
sequence of defect (hyper-) planes, characterized by very weak vertical bonds Ji ((ii) in Fig. 
2). Note that all physics of the second component, as well as all details of the microscopic 
implementation of the i interface between the subsequent layers, have been absorbed into two 
parameters: the directly measurable correlation length £ 2 and the nonuniversal amplitude 
Jo- 

A necessary condition for the reentrant dimensional crossover, described in the Intro- 
duction, is the smallness of the initial, bare value of the coupling in the temperature range 
around T c « T c \. 

Ji = J exp[-L 2 /6(T cl )] < k B T. (2) 

Under this condition, we can consider the coupling ([!]) as a weak perturbation with respect to 
the reference system consisting of noninteracting layers of the first component. Further, for 
the values of the RG parameter / such that ^2 C A' < Ii, typical fluctuations are correlated 
only over lengths much shorter than L 1; implying that the top and the bottom surfaces of 
each 1-layer are uncoupled. Provided the correlation length ^ of the first component has not 
been encountered yet, the system appears to be at the critical fixed point CFP S describing 
independent bulk critical behavior within the 1-layers coupled to the ordinary surface critical 
behavior at the noninteracting interfaces. Note that this scenario implies another condition, 
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Li»&(T cl ). (3) 

Both conditions, (0,0), require the layers to be thick. In the case of two different materials, 
T c i — T c2 = 0(1), so that £2(^1) = 0(02) and both inequalities are satisfied as soon as 
Li, L 2 ^> 02- The situation becomes less trivial in the case of a superlattice created by a 
weak periodic modulation of the properties of an originally uniform sample, then one has to 
be sure that the period of modulation is large enough to offset the smallness of the difference 
between the bulk Curie temperatures implying a relatively large £2(^1) oc (T& — T c2 )~ v , 
where v is the standard correlation length exponent in d dimensions. 

The crossover is now described as an RG-flow between three critical fixed points (Fig. 4): 
CFP S describing a sequence of uncorrelated thick layers of the first component ((ii) in Fig. 
2), with extensive properties dominated by the d-dimensional bulk critical behavior taking 
place inside the layers, while the associated ordinary surface scaling describes the observables 
localized at the surface. From CFP the system flows to CFP 1 describing the d' = (d — 1)- 
dimensional criticality in a system of uncorrelated thin films of the first component ((iii) 
in Fig. 2), and, finally, to CFP^ describing uniform <i-dimensional bulk behavior ((iv) in 
Fig. 2). The flow between the fixed points is driven by two scaling fields (cf. Fig. 4), 
the inverse thickness of the 1-layers, L~[ x , playing the role of a long wave length (infrared) 
cutoff in the ^-direction (cf. |16|), and the interlayer coupling strength Jj. Yet another 
one, scaled temperature field t± — (T — T cl )/T cl controls the departure from the critical 
manifold, ti c (Lx, Jj), containing the flow attracted by CFPb, towards the massive, high- 
and low-temperature fixed points. Instead of the scaled temperature field t\ characterizing 
the bulk criticality of the first component, one can use t = (T — T c )/T c defined relative to 
the observed critical temperature T c (Lx, J%) of the superlattice. As we will see below, the 
scaling forms look simpler when expressed via t, but lack any information about the shift in 
the critical temperature with respect to T c \. In the RG approach the scaling fields become 
functions of the logarithmic length scale /. While the explicit calculation of the crossover 
scaling functions requires application of non-perturbative methods such as the one used in 
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|TT|| , the leading singularities are determined by the RG flow in close vicinities of the fixed 
points. 

Specifically, at CFP S the linearized RG flow equations are 

dh/dl = (l/i/)ti, (4) 
d{l/L x )/dl = 1/L 1} (5) 
dJjdl = {d 1 - 2u!)Ji = {juMJi. (6) 

Here v is the standard correlation length exponent of the bulk <i-dimensional universality 
class. The second equation (|5]) reflects decrease in L\ expressed in units of the running 
inverse length scale A'; the equality of the RG eigenvalue in (||) to one reflects asymptotic 
isotropy of the bulk critical behavior described by CFP S . The RG eigenvalue for Jj is read 
from (1). Because the spins entering (1) are located at the surfaces, the scaling dimension of 
the spin density U\ determining the RG eigenvalue of Jj is that characteristic of the surface 
spin density. Correspondingly 7 n is the exponent characterizing the susceptibility of the 
surface spins to a perturbation by a surface field @. The crucial point now is that since 
the surface spins are correlated more weakly than the bulk ones, the surface susceptibility 
usually does not diverge (see the estimates of surface critical exponents in [||]): ju < 0. 
Consequently J, is irrelevant or marginal at the first encounter with CFP. Thus starting 
with a small coupling Jj, we are guarantied that it remains small until the system arrives 
at CFP'. Note that this picture implies a seemingly paradoxical prediction: a single layer 
of weak bonds cutting through a bulk sample effectively decouples the two halves, with two 
independent ordinary surface fluctuations developing at the two sides of it. 

The two other scaling fields, t\ and L^ 1 are relevant at CFP S . As Li{l) is positive 
definite, it is convenient to exclude I from the equations (^H), rewriting them as 

dln(ti)/dln(Li) = -I/ -1 , (7) 
d]n(Ji)/dln(Li) = -7UI/- 1 . (8) 

The crossover to the (i'-dimensional fixed point CFP' occurs when L\(l) decreases to the 
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order of the microscopic length scale a±. As a result of this first crossover (indicated by 
subscript D, for "direct," below) the scaling field 1/Li becomes irrelevant and disappears 
from the consideration, while the other two are renormalized to 

tin = h/r D (9) 
J lD = J l r D 111 . (10) 

Here the temperature rescaling factor 

t d « {al/Ltfl" (11) 

conveniently characterizes the first crossover. We will reduce the number of nonuniversal 
parameters by defining td not through a\, but rather through the critical amplitude X\ = 
0(ai) characterizing the divergence of the bulk correlation length £i = Xi(— t\)~ u at the 
low-temperature side of the bulk critical point T&: 

r D = {XyjUflr (12) 



In the RG formalism |T3[ the correlation length ^ characterizes the crossover from the critical 
to noncritical RG fixed points at nonzero values of t\. Therefore td determines the width 
of the temperature domain around T cl in which the dimensional crossover CFP S — > CFP' 
actually occurs: for t\ > td the ti-field becomes large and drives the system away from the 
critical manifold before 1/Li grows large; as a result the system never reaches CFP'. 

If, on the other hand, t\ ^ tjj the system arrives at CFP' with the values of the two 
relevant fields estimated by t\ ~ t\D, J% ~ JiD- In the absence of Jj (which is small 
at this stage) t\ completely determines the flow. The critical separatrix going into CFP' 
corresponds to a certain initial value tm = t' c = 0(1). One expects t' c < 0, as finite thickness 
suppresses ordering. In fact, due to our definition of 7~d via the bulk correlation length £i, 
the parameter t' c is universal: criticality in a free film of the first component of thickness Li 
happens at a temperature T' C {L\) at which the ratio takes a universal value 

Li/UK) = (-t' c T « 2.89, (13) 
11 



where the numerical estimate has been obtained by series expansion methods [|17]]. Expand- 
ing around the separatrix one obtains 

dfci/dZ = (l/i/)(fi-t / J ( 14 ) 
dJi/dl = {d! - 2u/)J< = {i/v')Ji, (15) 

where the prime marks exponents related to the (f -dimensional criticality. These equations 
are to be solved with the initial data given by t\ = tip, Ji = Jj£>. As a result of the first 
crossover the surface spins in (fl|) are strongly correlated with the rest of the corresponding 
layer, thus becoming <i'-dimensional in nature. Correspondingly, the scaling dimension of 
the spin density U\ in the RG flow equation for Ji is changed from the ordinary surface value 
uj\ to the cf -dimensional u>'. Correspondingly nonpositive 7x1 is changed to positive 7', so 
that Ji is a relevant perturbation at CFP' driving the second crossover back to CFP. The 
following analysis of the reverse crossover (as indicated by subscript R below) essentially 
repeats the one performed for CFP S — > CFP' above: we divide one of the two linear RG 
equations by another to obtain (the coupling strength Ji is positive definite) 

d\n(t 1 -t ac )/dlnJ i = l/j'. (16) 

A nonuniversal amplitude Jq = 0(k B TD~ d ') is defined so that at Jj(/) > Jq the d'- 
dimensional hyperplanes become strongly coupled and the system crosses over to the uniform 
bulk behavior at CFPb. As this happens, t\ is rescaled by the factor tr = (Jin/ 'Jq) 1 ^ , i. e. 

tx(l) = tin = (t 1D - t' c )/T R , (17) 

where (recall that 711 < 0) 

tr = (Jo/Jo) 1/7 '(^i/^i) H7lll/7V exp(-L 2 /7'6). (18) 

Just as previously, for ^> 1, i. e. \tio — t' c \ ^> tr, the system flows away from the 
critical separatrix before the reentrant crossover to CFP can take place. Hence tr measures 
in units of tio the width of the temperature domain in which the reentrant scaling can be 

12 



observed; in the original units of t\ the width is given then by the product t^tr. Within 
the domain t\R < 1 the system flows into the close vicinity of CFPf,. The critical separatrix 
is again defined by a certain (this time positive: interlayer bonds enhance ordering) value of 
t\R = t c R = 0(1). In fact, as at this point we have no other scale to measure the interlayer 
coupling amplitude Jo, we can fix the arbitrary constant Jq by setting t c R = 1. Having 
grown to the order of Jq the interlayer coupling amplitude Jj is absorbed into the rescaled 
value of the temperature field ti(l) = tm and becomes irrelevant. The reentrant CFPf, is 
a regular critical fixed point with temperature (in the absence of magnetic field) being the 
only relevant perturbation measuring the deviation from the critical separatrix (Fig. 4). 
The evolution of the latter is again described by the RG flow equation (£|) which has to be 
solved with the initial condition 

*i(0 = tm - 1 = (h/r D - Q/t r - 1. (19) 

Having established the principle features of the RG flow we are in a position now to de- 
velop a scaling description of the observable quantities. The two temperature scales tr> and 
t r defined by (2) and (3), together with the original value of the scaled temperature of the 
first component, t\, conveniently parametrize the scaling functions. If the bulk properties 
of the two components are known, then the only new nonuniversal parameter appearing in 
our description of the superlattice is the dimensionless ratio Jq/ Jq , essentially characterizing 
the strength of coupling between the spins across the interface between the two components 
(see the Appendix for more precise definition). The standard two-scale factor universality 
fj~8|| of the bulk critical points is thus extended to what may be called 2 + 2-factor univer- 
sality: all scaling functions of a superlattice are universal apart from the two independent 
bulk critical amplitudes of the first component plus the amplitude Jq/ Jq characterizing the 
interface and the bulk correlation length of the second component ^(Jci)- in fact, the two 
additional parameters enter the description via the initial value of the effective coupling 
Ji(l = 0) = Jq exp [7^2/^2(^01)]. Thus, if one were not interested in the (singular) dependence 
of the critical behavior of the superlattice on the thickness L2, only one extra amplitude, Jj, 
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would have to be added to complete the description of critical behavior in a superlattice. 
However, since for large L 2 the amplitude Jj is going to be anomalously small, and since L 2 
is easily measured, we prefer to split Jj into the nonsingular amplitude Jo and the singular 
exponential factor exp(— L 2 /£ 2 ), adding the bulk correlation length £2(^1) to the list of 
empirical parameters. 

We start by analyzing the shift in the critical temperature. Unfolding back the two 
renormalizations of the temperature field, one obtains from fl4]Jl9D 



T c = T cl [l + r D {t' c + t r )} = T 6 '(Lx) + T cl r D r R . (20) 

The shift consists of a larger, 0{L\ X ^ V \ shift to lower temperatures, slightly corrected by a 
much smaller, 0[exp(— Z^/^)], shift in the opposite direction. The last expression repre- 
sents the latter as a shift with respect to T' c (Li) = T cl [l + r D t' c ], the critical temperature in 
a film of the first component of thickness L\ with free boundaries. If this latter temperature 
is known, then the observed shift towards the higher temperatures, T C — T' C {L\), can be used 
to estimate the unknown amplitude Jo/ Jq entering the definition of r R (|TE|). 

Let us consider now an extensive observable, say the specific heat per unit volume C 
characterized by critical exponents a and a' in dimensions d and d! correspondingly. Since 
the surface-to-volume ratio in a superlattice vanishes as L\ t 2 — > 00, the observed signal at 
the first stage of the crossover is dominated by the interior of the 1-layers: as far as extensive 
properties are concerned CFP S is equivalent to CFPf, and the RG flow is indeed reentrant. 
This flow can be represented by the scaling form 

C{t x ) = (Li/ L)Ai_\ti\ ~ a Ax (ti D , tm) , (21) 



where tip, t\ R are given in (pjfi7|), and Ai_ is the bulk critical amplitude of the first com- 
ponent on the low-temperature side of the criticality. The factor {L\j L)Ai_ represents the 
amplitude of the contribution to C from the interior of the 1-layers at the first stage of the 
RG flow (governed by CFP S ) below T cl ; the choice of Ai_ i instead of A x+ seems natural in 
view of T c < T c i. The scaling function Ai is completely universal. It has the following limits: 
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it approaches 1 at tm — > — oo and the universal amplitude ratio Ai + /A\- at tm — > +00; 
at — >• it develops a singularity Ai oc I^idI" which cancels the one generated by the 
\t\~ a factor in (plj). At |tio — <C 1 one has At ~ — icl~ a '£i(^ifi)- The asymptotic 
behavior of the new scaling function C\ follows the same logic as the one employed above 
for A\. it takes finite limits at t\ R — > ±oo, develops a singularity, Ci oc at t\ R — > 

to compensate for the singular prefactor, and, finally, diverges as C\ oc |ti# — l| _a , when 

tin -> I- 

As already mentioned above, because in the Ising universality class the critical tempera- 
ture shifts satisfy TpT R <C T£>t' c C 1 a simpler scaling form is achieved in terms of the scaled 
temperature t = (T — T c )/T c — t\ — r D (t' c + r R ) shifted to the observed critical temperature 
T c (Li, Ji). Defining tp = t/ru, t R = tn/T R one can write 

C = (L 1 /L)A 1 _\t\- x A(t D ,t R ), (22) 

where the universal scaling function A takes the following limits: 

A « Ai+/Ai-, atl«i « t R , (23) 

.Awl, at t R < t D < -1, (24) 

A « C D± t~ a t a D - a \ at |t D | < 1 < |t fl |, (25) 

A « C R± r a t a D ~ a 't a R - a , at |t D | < |t fl | < 1. (26) 

This form has a simple graphical interpretation given by the double logarithmic plot of 
Fig. 3. The continuity of C{t) requires that both universal amplitude ratios Cd±, C r ± 
(remember that the scaling function has been normalized by the bulk amplitude A\J) are 
numbers of the order of unity. The last asymptotic expression in ( p6|) shows that while 
the ultimate divergence A = A±\t\~ a is characterized by the bulk exponent a, the critical 
amplitudes 

A ± = (L 1 /L)A 1± C e± t%-°' (27) 

are shifted on the logarithmic scale from the bulk amplitudes A\± of the first component by 
(a' — a) \n(l/r R ) (see Fig. 3), i.e. 
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HA±/A 1± ) = i_ [L2/6(Tc) + ( /v) ln(Li)] + (2g) 
a' — a 7 

More information about the critical state of a super lattice governed by CFPb can be ob- 
tained by studying the anisotropy of correlation functions. The basic ratio X = £^(T)/£||(T), 
of the critical amplitudes of correlation lengths across and along the layers, is given by expres- 
sions identical to (p7|), (p8|), where one has to use anisotropic correlation length exponents: 
u± — i/n — v at CFP S , CFPfy, but u± = 0, Un = v' at CFP', since at CFP' the correlations 
grow only in the direction along the layers. The result is: 

X~-4- (29) 

Similarly, at T c , the amplitude of the spin-spin correlation function, G±(z) oc z~ 2uj at z ^> L 
is suppressed compared to that of G\\ (x) oc x~ 2uj by an amplitude ratio 

n = GjG\\ ~X 2 " ~rjf. (30) 

One thus arrives at the large-scale description of the critical state of the superlattice as a 
highly anisotropic realization of the bulk universality class of the first component ((iv) in 



Fig. 2). In fact, following the standard ideas of scaling |19[ one can relate all anomalous 



changes in the critical amplitudes (0), (0) to the anisotropy X of the basic length scales. 

III. PLANAR ISING SUPERLATTICE: EXTENSION TO THE CASE OF 
ZERO-TEMPERATURE CRITIC ALITY IN D' = 1 DIMENSIONS AND 
COMPARISON TO EXACT RESULTS 

Despite the generality of the above consideration, of all the 0(n) spin lattice models in 
three and two dimensions only the d = 3, n = 1-superlattice behaves strictly according to 
this scenario, as the Ising model has regular finite-temperature critical fixed points without 
marginal operators in both d = 3 and d' = 2. The exactly solvable problem of a d = 2-Ising 
superlattice requires certain modifications because CFP' describing the one-dimensional 



Ising model is a zero-temperature fixed point. This results |TlJ| in a large shift in the 
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critical temperature, T c \ — T c = 0(1). More precisely, in the interesting case of a thick 
layer superlattice |TT| , criticality happens when the correlation length £ 1; characterizing 



the low-temperature phase of the first component, is much smaller than the width of the 
1-layers, L\. This means that the RG flow in this case comes to CFP' not from CFP S , 
but from the low-temperature fixed point LTFP describing the ordered phases of the first 
component. Nevertheless, many results can be formulated in terms of just one nonuniversal 
parameter characterizing LTFP: the linear energy Si of a domain wall separating spin-up 
and spin-down phases of the first component. 

Let us start our consideration from a close vicinity of T c \, where £1 ^> L\ and the 
crossover CFP S — > CFP' proceeds directly, without encountering LTFP. After leaving 
CFP S the scaled temperature ti is renormalized to t\o = h(Li/ 'X\-p v . Since v = 1 in 
the planar Ising universality class, below but close to T cl one can write t\z> = — (note 
that the normalization \t\\ = l/£i, making X\- = X\ + = 1 has been adopted in [|TT|j). 
Further, j n = in this universality class [|J, so that the coupling Jj oc exp(— is 
not renormalized near CFP S to the linear order in Jj. While one should expect logarithmic 
corrections to J, for this marginal case, those do not seem to affect the essential features of 



the exact solution. As discussed in |TTJ, at £i <C 1 



\t lD \=L 1 /^ l = 2L 1 E 1 /k B T = 2J'. (31) 

The last identity, in which £i(T) is the linear free energy of a domain wall between the two 
low-temperature phases of the first component, represents the dimensionless coarse-grained 
spin-spin coupling J' of the d = 1-Ising model onto which a strip of the first component is 
mapped after the direct, bulk-to-film crossover. At CFP' this spin-spin coupling transforms 
according to dJ'/dl = —1/2, or 

dtxjdl = 1 (32) 

with the critical separatrix corresponding to the extreme value t' c = — oo, as implied by 
T' c = in the d' = 1-dimensional Ising class. Note that the fugacity ( = exp(— 2 J'), taking 
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the initial value exp(ii£>), behaves at CFP' as a regular thermal scaling field with the RG 
eigenvalue equal to one. Consequently, the criterion for CFP'-dominated behavior is C ^ 1, 
or, equivalently, tm <C — 1. 

Meanwhile, magnetization scales trivially, u' = 0, at CFP'. Thus the eigenvalue of 
Jj is d! — 2a/ = d = 1 and the corresponding RG equation at CFP' can be rewritten as 
dlnJi/dl = 1. Dividing one RG equation by another we obtain 

dt 1 /dkiJ i = l, (33) 

from which one can construct a renormalized temperature 

t 1R = t 1D - ln^/J*) = L 2 /£ 2 - Li/6 + 0(1). (34) 

It is convenient to define Jq here in such way that the criticality occurs at tm = 0. Then, 
in full agreement with [11], up to corrections vanishing as 0(Li \), the critical temperature 



T c of the superlattice can be obtained from the equation 

L 1 /^(T e ) = L 2 /UTc)=9c i (35) 

where the last identity defines what may be called the scaled thickness of the layers. 

If this equation is satisfied at g c <C 1, the case called the thin layer limit in ||11|| , then the 



flow indeed never encounters LTFP. However, in this case the initial value of the coupling 
Ji oc e~ 9c is not small, so that the criticality is reached before tm = —g c becomes large and 
negative so as to display any features of the d! = 1 critical behavior. The flow in fact never 



leaves CFP. The reentrant crossover does happen in the opposite, thick layer limit []TT 
g c 3> 1, however in this case the behavior of the system on length scales between £i and L\ 
is governed by the low-temperature, non-critical fixed point LTFP. 

To describe criticality in this limit we start at scales larger than £i(T) well below T c . 
If the RG is implemented with rescaling the block spin by the factor of the area of the 
block (rather than the square root of the area leading to the Gaussian fixed point described 
by the Ornstein-Zernike functional |20|), then at these scales fluctuations in the amplitude 



of magnetization are completely suppressed. The only allowed fluctuations are domain 
walls separating domains of two different spin orientations. At these scales the walls look 
geometrically sharp and are characterized (neglecting lattice effects) by a single parameter: 



free energy per unit length Ex. As discussed in [EI], a single strip of the first component 
becomes identical to a d = 1-Ising model: spins are fully correlated across the strip; a 
domain wall cutting across the strip has finite free energy E1L1 playing a role of an effective 
ferromagnetic coupling between neighboring blocks of size £1 x L±. Characterization of the 
effective d = 1-model is completed by specifying the value of a block spin, mi(T)Li£i, where 
m\(T) is the bulk magnetization density in the first component. 

The crossover from CFP' to CFP^ proceeds exactly as described in the beginning of 
this section, because the representation fl3T|) of tm via the linear surface energy Ei is valid 
independently of whether the system has arrived at CFP' from CFP S or from LTFP. After 
the reverse crossover the thermal field is renormalized to 

t 1R = t 1D - ln^/J*) = L 2 /6 - L^/ksT + 0(1), (36) 

leading to the universal relation for the critical temperature 

Si(r c )fr(:r c ) = u 

holding up to corrections of order 0(a/L). Note that the relation Ei = kgT '/2£i used in ([31]) 
is valid in the planar Ising universality class up to corrections of order 0(a/£i), which were 



neglected in the field-theoretical description of Ref. [fnj| . The expression (|37|) however has 
a much wider range of validity: it requires only that the strips are thick, with no condition 
on the values of the correlation lengths in the two components. 

For an extensive observable like the specific heat per unit area we can write now 

C = ~A'(T) exp(-a't 1D )A(t 1R ). (38) 

Here the exponent a', equal to —1 in the case of the specific heat, and the critical ampli- 
tude A'(T) characterize singularity C = A'(~ a ' in the one- dimensional regime. The scaling 

19 



function A has the following limits: it approaches 1 at tm ^> 1; diverges as A ~ C±t\R at 
txjj "C 1 with the amplitude ratios C± = 0(1) being universal; an exponential singularity 
A oc exp(a'ti^) cancels the singular prefactor in (|3"8] ) at £nj — > — oo. The critical divergence 
of the specific heat, 

C « ~A'(T C ) exp(-a^ c )C±trfl = (39) 
is characterized by the amplitude 

A~ A'(r c )exp(-a^ c ), (40) 



where we have taken into account that dtm/dt = O(L) cancels the L _1 prefactor in fl3"9|). 

In order to complete the estimate of A we have to know the amplitude A'(T C ). The one- 
dimensional amplitudes of extensive observables (per unit length of a layer) can be generated 
from the well-known expression for the free energy density /' of a one- dimensional chain at 
temperature T and magnetic field H (see and references there) 



f/k B T ~ ((HmM 2 + (C/^i) 2 ) - , (41) 

where we have taken into account that £i plays the role of the effective cutoff for the one- 
dimensional behavior along the layers, and both £i and mi are taken at T = T c . Differenti- 
ating with respect to the temperature T we obtain (at T = T c and H = 0) A'(T) = g^B/Ci- 
Combined with (fiol) this gives 

A ~ fc B ^e- fl ^r 2 . (42) 



This agrees fully with the result of [11]]. Other amplitudes can be similarly estimated and 
agree with the exact solution. 



An interesting feature of the form (38) is that below T c , at tm <C — 1 the exponential 
part of the scaling function A compensating for the (now spurious) singular prefactor, leads 
to 

C oc exp(d / L 2 /e 2 ). (43) 
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This can now be recognized as pointing to the dominance of the fluctuations in the second 
component as T c2 is approached: in fact, due to the self duality of the planar Ising model, 
a completely similar construction can be developed starting from low-temperatures and 
focusing on the dual variables in the 2-layers. Both scenarios lead to essentially the same 
results. 



IV. HEIZENBERG AND XY SUPERLATTICES 

The other 0(n) models can be analyzed similarly; only the most interesting results will 
be briefly discussed below. 

The case of the three-dimensional Heizenberg (n = 3) model is very close to the planar 
Ising one, as the (d' = 2, n — 3)-class is believed to be characterized by a zero-temperature 
critical point with the effective coupling being marginally relevant f2~2f . This effective cou- 
pling in a single Heizenberg film of the first component is the dimensionless helicity modulus 
p3| T[. In the low-temperature phase of the d = 3-Heizenberg model the helicity modulus 
T 1 has the dimension of inverse length, effectively defining the correlation length £1 = T^ 1 . 
In the thick layer limit the critical temperature T c of the superlattice is again shifted well 
below T c \ of the first component. Thus on scales smaller than L\ the first component is well 
ordered and the additive approximation V = YL\ ^> 1 is well justified [p3[ . Here F± can 



be (at least in principle) measured independently in a bulk sample of the first component. 
Once the direct crossover to CFP' has occurred, one can form = k^T ' /V and use the 



well known d = 2-Heizenberg model RG equation p2 

dt 1D /dl = l/2vr, (44) 

while the coupling Jj is trivially additive 

dJ./dl = d'Ji, (45) 

with d! = 2 here. Combining these two equations we obtain, in complete analogy with ( P7| ) , 
a universal expression satisfied at the critical point of a Heizenberg superlattice: 
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k B T c AttLi 

Quite analogously to the consideration given above for the planar Ising case one arrives at 
the picture of the critical temperature shift T c \ — T c = 0{T C \ — T c2 ), large scale fluctuations 
between in the wide temperature range between T cl and T c2 , and exponential dependence 
of critical amplitudes on the thickness of the layers. 

The case of a XY (n = 2) superlattice in d = 3 is special because, on one hand, 
the d! = 2-class is characterized by a finite-temperature Kosterlitz-Thouless transition, so 
that we expect the shift in the critical temperature to vanish in the limit of (atomically) 
thick layers. On the other hand, the spontaneous magnetization is zero in d' — 2 at any 
finite temperature. Thus the coupling Jj through the 2-layers is crucial in supporting long- 
range order in a three-dimensional superlattice and we expect the observed spontaneous 
magnetization to be anomalously suppressed when Jj oc exp(— L 2 /£ 2 ) is small, even well 
below T c pd T cl . The suppression is easily estimated by combining the RG equation for the 



order parameter dm/dl = —lu'tti, where the scaling dimension |22| 



J = ksT/An^ = ksT/ATrT^, (47) 



with the additive growth of Jj given by fl45|) ; note that in the case of a superfluid the helicity 
modulus T is proportional to the superfluid density p s . The result is 

mi(T) 87rri(T)£ 2 (T) L x 

Note that k B T /T\ essentially plays the role of £i in this universality class, as well. One 
can thus see that in the case of L1/L2 = 0(1), spontaneous magnetization of a superlattice 
decreases exponentially as soon as £1 exceeds £ 2) l° n g before £1 becomes comparable to the 
thickness L\. It would be interesting to see if this behavior could be observed at a superfluid 
transition of liquid helium-4 in some matrix with periodically modulated properties. 

Similarly, a d = 2 — XY superlattice may be realized in a helium film deposited on a pe- 
riodically modulated substrate (the latter could be itself cut out of a solid three-dimensional 
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superlattice). In this case the d! — 1 Xl^-model orders only at T = 0. Thus we expect T c of 
the superlattice to be shifted well below T cl . At such temperatures the 1-layers are equiv- 
alent to a sequence of one-dimensional (classical) AF-chains, described by the Gaussian 
Hamiltonians 

Uj = J dx-F^Vfa) 2 , (49) 

where is the phase of the XY order parameter, j labels the layers. The effective one- 
dimensional helicity modulus T' is given by the product TxLi, as in all other cases considered 
in this section. Each pair of neighboring layers is coupled through the 2-layers separating 
them via 

Ttjj+i = Ji J dx cos(0j —(f)j + i). (50) 

Combining the flow of the one-dimensional Gaussian coupling dV jdl = —V with the ubiq- 
uitous (pE5]) (ef = 1 here) we arrive at the condition of criticality YiLiJi/ \k B T) 2 = 0(1), 
which with logarithmic accuracy can be rewritten as 

L 2 /UT c ) = HLi) + 0(l). (51) 

A little thinking shows that for L\/ L 2 = 0(1) the critical temperature of the superlattice is 
in fact shifted closer to the lower of the critical temperatures of the two components; it stays 
however sufficiently above T c2 so that £2 remains much smaller than L 2 and our focusing on 
fluctuations in the 1-layers remains well justified. 

V. SUMMARY 

In summary, this paper addresses the following, as yet hypothetical, experimental project: 
A sequence of superlattices is built of the same two ferromagnetic components, but with 
different basic thicknesses L\ and L 2 (Fig. 1). The bulk critical properties of the components 
are assumed to be well known, the thicknesses Li,L 2 3> a, i.e. large compared to any 
important microscopic length of the two materials, and no long-range dipolar interaction 
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is present. Given the well developed theory of critical phenomena in the Ising, XY, and 
Heizenberg universality classes, can one predict the variation of the critical properties of the 
superlattices with the changes in L\ and L2? 

Naively one may expect that in the limit of L%, L 2 — > 00 all extensive properties of the 
superlattice can be expressed as a weighted average of the bulk densities of the components. 
Indeed, if the limit of infinite thicknesses is taken at fixed temperature T 7^ T c , surface 
corrections to the bulk density average vanish as (£1 + £2)/-^, where the correlation lengths 
£1, £2 P^y the role of penetration depths for the perturbations induced by an interface in 
the two adjacent layers. However, experiment on a given sample studies the T — > T c limit 
at fixed thicknesses L\, L 2 , so that the correlation length(s) are bound to exceed the period 
of the superlattice making the surface contributions crucial. The thick layer limit for the 
critical properties of a superlattice is singular, the thicknesses Li,L 2 have to be combined 
with certain functions of the temperature to form proper scaling fields in terms of which a 
consistent scaling description is achieved. 

Such scaling description has been constructed above based on the qualitative picture of 
reentrant dimensional crossover [IT]. In this picture (Fig. 2) a superlattice, at its critical 



temperature T c , displays three different asymptotic types of critical behavior, when probed 
on different length scales. These three asymptotic models are, in the order of increasing 
wave length, (ii) a sequence of noninteracting semi-bulk samples, exhibiting ordinary surface 
critical behavior at the surfaces; (iii) a sequence of noninteracting thin layers, exhibiting 
d! = d— 1-dimensional behavior; (iv) a highly anisotropic, but uniform critical system of the 
same universality class as that of the components. Away from T c the finite correlation length 
falls into one of the three regimes, so that instead of probing different length scales at T c 
one can study the crossover by observing the asymptotic behavior at different temperatures 
approaching T c . 

Each of these limiting types of behavior corresponds to a previously studied renormaliza- 
tion group fixed point. Combining those into a global RG flow (Fig. 4) we have constructed a 
rather complete scaling picture of critical phenomena in superlattices composed of thick lay- 
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ers of two components. The simplest case is that of a d = 3-Ising superlattice, characterized 
by the possibility of a finite-temperature phase transition in a single layer of the first compo- 
nent. As a result, in the limit of thick layers the critical temperature T c of the superlattice 
approaches that of the stronger of the two ferromagnets, T c \. The two crossovers then hap- 
pen in a (relatively) narrow vicinity of T c . They are conveniently characterized in terms of 
the reduced temperature scales td, tr given by (0), (|18|)- These scales represent (Fig. 3) the 
relative widths of the temperature domains in which the superlattice behaves, respectively, 
as models (ii) and (iii) above. The simplest scaling description ( p6|) is then achieved in terms 
of the reduced temperatures tp = (T — T c )/te>T c and tR — (T — T c )/trTdT c , shifted to the 
observed critical temperature T c . As represented graphically in the double logarithmic plot 
of Fig. 3, the ultimate divergence is characterized by the critical exponents characteristic of 
the bulk, d-dimensional behavior. However the presence of the intermediate domain of the 
cf-dimensional behavior (model (ii)) leads to anomalous changes in the critical amplitudes, 



as seen in (Fig. 3) and expressed analytically through ( [H] , [[SD in (|28|). 

The situation is somewhat more interesting in the other three- and two-dimensional 
spin-vector universality classes, where the weak bonds between subsequent 1-layers play a 
crucial role in maintaining the magnetic long-range order in the system. This leads to a large 
shift in the critical temperature of a superlattice, as implicitly given for various classes by 



the universal expressions (|37|j46| , fl9|), as well as the anomalous, exponential suppression in 
the amplitude of spontaneous magnetization (f48|) in a three-dimensional XF-superlattice. 
Where applicable, the results obtained here by the method of this paper agree with the 



analysis of the exact solution for the two-dimensional Ising superlattice |T(| given in |TT 



The three-peaked form of the temperature dependence of the specific heat, characteristic 
of that model generalizes to the three-dimensional Heizenberg class. It has a transparent 
physical interpretation: most of the spin degrees of freedom in such a superlattice get ordered 
at one of the bulk critical temperatures of the two components, leading to finite, but large- 
amplitude peaks at both. However, a small, 0(L _1 ), fraction of the degrees of freedom 
participate in strong large-scale fluctuations at all T c2 ^ T < T&, giving small- amplitude, 
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but diverging contribution to the specific heat at T c . 

The abundance of interfaces in a three dimensional magnetic superlattice may greatly 
improve chances of experimental observation of surface and two-dimensional scaling. In par- 
ticular, the elusive surface scaling exponents enter the definition of the parameter tr in ([18]) 
and can be obtained, for instance, by analyzing the dependence of the critical amplitudes, 
(p8j) and Fig. 3, on the thickness of 1-layers L\. 

Finally, the main restriction of the present analysis seems to be the neglect of dipolar 
interactions. When present, these interactions will lead to demagnetization effects making 
Heizenberg universality class asymptotically unstable. Otherwise, most of the results pre- 
sented above can be easily modified in the case when the interlayer interaction dependence 
on the thickness of the separating layers, Jj(L 2 ), is a power-law rather than an exponential 
as was always assumed above. However, the assumption of local interactions enters our con- 
sideration in much deeper ways than just through the exponential form of Jj(L2), so that 
the required revision could be much more serious. Aniz-ferromagnets, of course, are free 
from the complications arising from the dipolar interactions. However, pinning of the spin 
density waves by the interfaces between the layers in a superlattice gives rise to a different 
source of the thickness dependent effects beyond the scope of this paper. In any case, the 
above consideration seems to represent a reasonable first step towards a complete analysis 
of scaling in superlattices. Note that recently critical phenomena in superlattices have been 
analyzed in the frames of the mean- field Landau- Ginzburg approximation ||24|| . A direct 



comparison of our results to theirs is hardly feasible in view of the failure of the mean-filed 
approximation d = 2 and 3 considered here. However a combination of the results of Ref. |24j 



with a renormalization group expansion in d = 4 — e and the scaling framework developed 
above looks promising, although computationally very hard, direction for future work. 
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APPENDIX A: INTERLAYER COUPLING 



Here we derive the effective spin-spin coupling (|T]) from the Ornstein-Zernike spin-density 



functional 





dz d d 'x'{s 2 + (6) 2 (Vs) 2 } 



(Al) 



2X2 Jo 



describing a layer of the second component on the length scales exceeding the correlation 
length £ 2 ; the other parameter entering ( |A1| ) is the bulk spin susceptibility \ 2 . This func- 
tional corresponds to the Gaussian RG fixed point describing the bulk paramagnetic phase 
of the second magnet. In a layer of finite thickness L 2 , the bulk part of the functional has to 
be complemented by the surface contributions, which at the paramagnetic fixed point take 
the form 



In the above equation z\ = 0, z 2 = L 2} the phenomenological parameter b accounts for the 
change in the number of nearest neighbors for the surface spins, while the surface fields hj 
are induced by the neighboring spins in the layers of the first component: 



where s 1 = and s 2 = s,f +1 , as defined after ([[]), while J 12 is the effective coupling 
across the interface. The minimization of the total free energy + F sl + F s2 proceeds 
via Fourier expansion in the x'-plane. For each Fourier harmonic of the spin density one 
obtains an independent one-dimensional functional. Minimizing the latter and summing all 
contributions one obtains in the limit exp(— L 2 /^ 2 ) C 1 an effective interaction between the 
surface spins of the first component 




(A2) 




(A3) 




(A4) 



k' 




(A5) 
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Since the Fourier components of the spin density with wave numbers k' > £f has been 
already eliminated by the RG transformation, the factors (1 + ~ 1- Using ( |A3|) 

one finally obtains the interaction of the form (0) with the amplitude 



26/X2 
[6 + 6/X2] 



J o = r L , ,» 12 J 12- ( A6 ) 
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FIGURES 



FIG. 1: Schematic view of a superlattice in, a), d — 3 and, b), d — 2 spatial dimensions. 

FIG. 2: Graphical representation of the reentrant dimensional crossover. The observed 
behavior of the system changes from a noncritical superlattice, (i), to the sequence of thick 
layers of the first compnent, (ii), connected by exponentially weak links mediated by non- 
critical fluctuations in the second component. From this limit, the direct crossover takes the 
system to a sequence of weakly coupled thin layers, (iii), while the reentrant crossover takes 
it back to the bulk, although highly anisotropic, critical behavior, (iv). 

FIG. 3: Schematic dependence of the specific heat per unit volume C on the scaled 
temperature t = (T — T c )/T c in a three-dimensional Ising superlattice. 

FIG. 4: Reentrant renormalization group flow in the critical manifold of the parameter 
space representing the d = 3 Ising superlattices. The three fixed points, defined in the 
text, are connected by critical separatrices (bold lines). They correspond, in the order in 
which they are encountered by the flow, to the limits depicted by (ii), (iii), and (iv) in 
Fig. 2. The direct crossover, from CFP S to CFP', and the reentrant one, from CFP' to 
CFPf,, are driven by two scaling fields, measured by the inverse thickness 1 and by the 
interlayer coupling strength Jj, respectively. Note that yet another relevant scaling field, t, 
directed roughly perpendicular to the plane of the figure, drives the system away from this 
two-dimensional critical manifold to the high- and low-tmeperature fixed points. 
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